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ABSTRACT 

Context. The Great Observatories Origins Deep Survey (GOODS) has provided us with one of the deepest multi-wavelength views of 
the distant universe. The combination of multi-band photometry and optical spectroscopy has resulted in the identification of sources 
whose redshifts extend to values in excess of six. Amongst these distant sources are Lya- emitters whose nature must be deduced by 
clearly identifying the different components that contribute to the measured SED. 

Aims. From a sample of Lya emitters in the GOODS-S field with uncontaminated photometry and optical (red) spectroscopy, we 
select a spatially compact object at a redshift of 5.563 (Lya) that shows a second emission line, identified as Niv] 1486 A. The SED 
is modelled in a way that accounts for both the N iv] line emission and the photometry in a self-consistent way. 
Methods. The photoionization code CLOUDY is used to calculate a range of nebular models as a function of stellar ionizing source 
temperature, ionization parameter, density and nebular metallicity. We compare the theoretical and observed magnitudes and search 
for the model parameters that also reproduce the observed N iv] luminosity and equivalent width. 

Results. A nebular model with a hot blackbody ionizing source of around 100 kK and a nebular metallicity of ~5% of solar is able 
to fit the observed SED and, in particular, explain the large apparent Balmer break which is inferred from the pure stellar population 
model fitting conventionally applied to multi-band photometric observations. In our model, an apparent spectral break is produced by 
strong [Om] 4959, 5007 A emission falling in one of the IR bands (IRAC1 in this case). A lower limit on the total baryonic mass of 
a model of this type is 3.2 x 10 s M Q . 

Conclusions. It is argued that objects with Lyo- emission at high redshift that show an apparent Balmer break may have their SED 
dominated by nebular emission and so could possibly be identified with very young starbursting galaxies rather than massive evolved 
stellar populations. Detailed studies of these emission nebula? with large telescopes will provide a unique insight into very early 
chemical evolution. 

Key words, cosmology: early Universe - galaxies: formation - galaxies: photometry - galaxies: starburst - galaxies: stellar content 
- ISM: abundances 



1. Introduction 

The Lya emission line is a beacon that can be recognised from 
the ground in objects with redshifts larger than about 1.8. When 
seen on its own, its interpretation is complicated by radiative 
transfer effects both within the emission region itself and in the 
intervening Lyman forest. With the detection of other emission 
lines and/or a recognisable continuum, a spectrum can how- 
ever provide a rich source of information about physical con- 
ditions and chemical composition. Of particular interest for the 
identification of primordial starforming objects are the so-called 
'dual-emitters': sources that emit both Lya and the Hen line at 



* Based on observations made at the European Southern Observatory, 
Paranal, Chile (ESO programme 170.A-0788) The Great Observatories 
Origins Deep Survey: ESO Public Observations of the SIRTF 
Legacy/HST Treasury/Chandra Deep Field South.); on observations ob- 
tained with the NASA/ESA Hubble Space Telescope obtained at the 
Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc.; and on obser- 
vations made with the Spitzer Space Telescope, which is operated by 
the Jet Propulsion Laboratory, California Institute of Technology under 
a contract with NASA. 



1640 A ( |Dawson et aL]|2004[ |Nagao et aL]|2005| |Ouchi et al" 
2008[ ) which are interpreted as very metal-poor nebulae ionized 



by massive, hot and presumably metal-free stars. So far, interme 
diate band imaging surveys have not been very successful in dis- 
covering genuine, stellar-ionized H/Hen dual emitters (Nagao 
|et al.|2008| >. 

The Great Observatories Origins Deep Survey (GOODS) 
spectroscopic observations in the Chandra Deep Field South 
( [Vanzella et aLpOOB] [2006] [2008 20091} [Popesso et al.|[2009 



have accumulated, using the FORS2 and VIMOS instruments on 
ESO's Very Large Telescope, spectroscopy of essentially all ac- 
cessible sources. The source GDS J033218. 92-275302. 7 (here- 
after GDS18.92-02.7) was selected for the spectroscopic pro- 
gramme as a V 6 o6-band dropout and has a blue /775-Z850 colour 
due to a strong Lya emission line falling in the /77s -band. It was 
subsequently studied in the context of SED fitting by Stark et al. 
((2007) and by |Wiklind et al.| ( [2008] ) who both conclude, on the 
basis of its apparently large Balmer break, that it is a high mass, 
old stellar population. Since the restframe UV images obtained 
with the ACS on HST show a very small angular size, the im- 
plication is that this is an example of an ultradense spheroid at 
a very early epoch. A significant fraction of other Lya emitters 
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from GOODS, showing apparent Balmer breaks, have also been 
interpreted as evolved, massive stellar populations (Pentericci 
et al.|2009) . 

The presence of strong emission lines in the UV spectrum 
carries implications about the effect of nebular emission pro- 
cesses on the rest of the observable SED. There are circum- 
stances in which the combination of nebular continuum from 
bound-free and two-photon emission together with emission 
lines can influence the appearance of, or even completely domi- 
nate, an underlying stellar population throughout the UV, optical 
and NIR spectrum. In this paper, we exploit the presence of a 
second emission line, identified as the [Niv], Niv] intercombi- 
nation doublet at 1483, 1486 A, to estimate the implied SED and 
compare with the extensive GOODS photometry of this object. 

The presence of significant N iv] emission is rather rare in 
known astronomical objects. At low redshifts it is seen in some 
compact planetary nebulae and in peculiar very hot stars such as 
j] Carinse (Davidson et aT7| | 1986| >. It is also sometimes detected 
in AGN and quasar spectra but with an intensity that is low with 
respect to lines such as C iv and N v, although there are rare ex- 
amples of 'Nitrogen Loud' quasars that exhibit unusually strong 
nitrogen li nes, e.g. [Baldwin et al.|p003) . At a redshift of 3.4, the 
Lynx arc (Fos bury et al.|2003| hereafter F03), is known to show 
N iv] with a strength comparable to other UV intercombination 
lines such as O in] and C m] and the resonance doublet C iv. This 
object is interpreted as an Hn region ionized by very hot stars 
(T e ff ~ 80 kK) with a nebular metallicity a few per cent of solar. 
In this case, the entire observed SED is nebular in origin and the 
stellar emission is not detected longward of Lya. Glikman et al. 
(2007) report two Niv] emitting low luminosity 'quasars' one 
of which, DLS 1053-0528 at z = 4.02, has an intensity ratio of 
N iv]/ C iv = 1 .9, a weak continuum and an SED which is con- 
siderably bluer than a typical quasar. This object has relatively 
narrow lines and may be more closely related to the Lynx arc 
than to an AGN. 

In this paper we examine a sample of Lya emitters from the 
GOODS spectroscopic data with redshifts spanning the range 
3.5 < z < 5.9 and present new photometric measurements from 
the most recent releases of the HST, VLT and Spitzer imag- 
ing data. For our selected Niv] emitter, we use a photoion- 
ization modelling code to investigate the range of model neb- 
ula? that will produce the observed N iv] emission flux. Given 
the GOODS broad-band photometric measurements that extend 
from the U-band to the Spitzer MIPS 24 fim channel, we cal- 
culate the contribution of the nebular line and continuum emis- 
sion to these bands in order to constrain the nature of any non- 
nebular light that may contribute. In addition to modelling the 
SED of this one object, we address the possibility that nebu- 
lar emission may contribute significantly to other sample mem- 
bers, a possibility that has already been suggested by Zackrisson 
et al. ( 2008| l. Our conclusions are of relevance to the interpre- 
tation of the SEDs of other high redshift Lyo- emitters, where 
the emission line strength suggests that the nebula may be a 
significant contributor to the spectrum longward of the Lya 
line. We adopt a cosmology with £2a = 0.7, Om = 0.3 and 
Ho = 70 km s _1 Mpc~'. The AB magnitude system is used 
throughout. 

2. Data 

2.1. Sample 

We have spectroscopically selected a sample of Lya emitters 
from the publicly available data for the southern field of the 
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Fig.l. Image cutouts of GDS18.92-02.7 from HST/ACS (top- 
row) and Spitzer/IRAC (bottom-row). The ACS stamps are 
5 arcsec on a side and the IRAC stamps are 25.6 arcsec. 



Great Observatories Origin Deep Survey (see |Dic kinson & 
GOODS Legacy Team 2 001 1 Giavalisco et al.|2004| for a review 
of the GOODS project). The multi- wavelength observations pro- 
vide us with information on the U (VLT), B435 V606 '775 and z&so 
(HST), Js,H,Ks (VLT), 3.6, 4.5, 5.8, 8.0 and 24 fxm (Spitzer) 
bands. Many other observations have been carried out, including 
deep integrations in the X-ray and radio domains from Chandra 
and the Very Large Array, respectively. 

The spectroscopic observations carried out with the ESO 
VLT FORS2 instrument yield a wavelength coverage from ap- 
proximately 0.55-1 yum with a resolving power of R « 660. This 
programme has produced about one thousand redshift determi- 
nations between redshift 0.5 and 6.2, with more than one hun- 
dred Lyman break galaxies confirmed at redshifts beyond 3.5 
dVanzella et al.|2005j [20061 [20081 . 

From the GOODS data product archive 
http://www.eso.org/science/goods we used the FORS2 
online engine to select all the qf (quality flag) = A 'emission' 
objects at z > 3.5. A few objects were excluded because of 
spectral contamination, one because it showed no Lya emission 
and one because of positional uncertainty. This process yielded 
23 objects from which an additional five were excluded due to 
large photometric errors in some passbands. 

2.2. Photometry 

Due to the different telescopes in use, the GOODS images are 
characterised by a large range in point spread function (PSF) 
dimension. An unbiased estimate of the colours is essential for 
analysis of the SED. We have used a PSF-matched photometric 
technique to obtain the AB magnitudes of each object in each 
band in a given aperture size. In the presence of potential con- 
tamination by neighbouring objects, we have used an optimum 
circular aperture which maximises the S/N of each object in all 
passbands simultaneously. For example, in the case of the match 
between ACS and groundbased Ks or IRAC images, compo- 
nents may be clearly resolved by ACS but be merged in IRAC, or 
very red objects may be detected at the longer wavelengths but 
show no counterpart in the optical images. A curve-of-growth 
analysis yielded the optimum uncontaminated aperture. 

Each source was examined in a 12 x 12 arcsec 2 cutout and 
a Monte-Carlo method used to estimate the background. The 



photometric measurements were made with Sextractor (Bertin & 
[Arnouts|T9 96 ), where the background used was consistent with 
Monte-Carlo simulations. Finally, by selecting some bright but 
unsaturated stars near each galaxy, we computed the convolu- 
tion kernel required to match the PSF of all the images to that 
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with the largest. We generated about 200 PSF kernels of stars in 
different filters in order to find the best kernel using a y 2 min- 
imisation method. All magnitudes were corrected to a 6 arcsec 
circular aperture. The results for all bands appear in Table|4] We 
also use Sextractor to measure the 'half-light' (or 'effective') ra- 
dius of each galaxy in the z^so-band (also listed in Table [4]). 

2.3. The Niv] emitter 

Of the 18 Lyo- emitters in our sample, only one shows a clearly 
detected second emission line. Galaxy GDS18.92-02.7 with 
Z850=24.5, was selected as a V606-band dropout and has been 
confirmed to be at redshift 5.563 (redshift of the Lya line uncor- 
rected for absorption). The spectrum, included in Figure [3] and 
shown in detail in Figure |4] was extracted from several individ- 
ual FORS2 multi-object masks and corresponds to a total expo- 
sure time of 14.4 ksec (Vanz ella et al.|20 06). Detailed measure- 
ments of the spectrum of this object are presented in Vanzella 
et al.| ( |2009b[ ) and the parameters relevant to our modelling are 
summarised in Table [T] 
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Fig. 2. FORS2 spectra and photometric SED - part 1 (redshifts 
3.6—4.8). The photometric points are: #435 V^ob '775 and z&so 
(HST), Js, H, Ks (VLT), 3.6, 4.5, 5.8, 8.0 (Spitzer). The sources 
are presented in order of increasing redshift and so can be cross- 
referenced to Table [4] The (red) + 'plusses' for each SED repre- 
sent the wavelengths of Lya, the 4000 A B aimer break and the 
[0 111] 5007 A line. The photometric error bars are lcr. 

The main spectral features are the Lya emission line (rest- 
frame equivalent width (EW) ~60 A), the break of the continuum 
just blueward of the line, and the intercombination emission dou- 
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Fig. 3. FORS2 spectra and photometric SED - part 2 (redshifts 
4.8-5.9). 
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Fig.4.FORS2 spectrum of source GDS18.92-02.7 showing the 
Lyor and Niv] emission lines and the ACS and VIMOS filter 
passbands. The expected positions of the N v, O iv and Si iv lines 
are also marked. 



blet [Niv], Niv] 1483.3, 1486.5 A, whose intensity ratio can be 
used to estimate electron density (see Figure|4]i. There is no clear 
detection of Si iv 1394, 1403 A in absorption or emiss ion nor of 
N v 1239, 1243 A in emission. IVanzella et all (|2009b) also note 



that the object is extremely compact but marginally resolved in 
the ACS images. 
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Table 1. Summary of the physical quantities derived from the 
optical spectral features of GDS18.92-02.7 taken from Vanzella 
|etaLld2069bb. 



Quantity 



Value 



Comment 



z ( Lya peak) 
L( Lya ) 
EW ( Lya ) 
z(Niv]) 
L(Niv]) 
EW (Niv]) 
L(Sirv) 
L(Nv) 



5.563 

3.8±0.3xl0 43 



erg sec 



59™ A(restframe) 
5.553 (1486.5 A) 
1.3+^xl0 43 ergsec-' 
22+% A (restframe) 

< 0.11 ± 0.07 x 10 43 erg sec 

< 0.11 ± 0.06 x 10 43 erg sec 



no abs. correction 

89 A with Ctnipho, 
high-density limit 

33 A with Ctniphot 
1 (1393.8, 1402.8 A) 
1 (1238.8, 1242.8 A) 



Equivalent width values are also given using the continuum level 
derived from the photometry (Ctm phot ). 



3. SED modelling 

In this section, we first discuss the modelling of 
GDS 18.92-02.7. We then use a grid of similar photion- 
ization models in order to assess the potential importance of 
nebular emission for interpreting the SED of the other sources 
in our sample. 



3. 1 . GDS J03321 8. 92-275302. 7 

Two previous studies have attempted to fit the SED of 
GDS 18.92-02.7 using the apparently large Balmer break to re- 



solve the age-extinction degeneracy. Stark et aL] ( |2007| , using 
ID=32_8020, and |Wiklind eTaT1 ( |2008| l, as BBG 5197, both ob- 
tain an age of 0.9 Gyr, a stellar mass of close to 10 n M o and 
zero extinction. The latter authors, who use the two longer wave- 
length IRAC channels, 5.7 and 8.0 fim, note that their fit is wors- 
ened by a deviant 5.7 ftm point. 

Our nebular modelling is motivated by the realisation that 
a source with such strong emission lines in the FUV will nec- 
essarily have nebular contributions at longer wavelengths. Also, 
by analogy with the Lynx arc at the lower redshift of 3.4, it is 
possible that the significant jump in brightness between the Ri- 
band and the IRAC 3.6 /im filter could be caused by strong [O m] 
4959, 5007 A emission lines producing an apparent, but false 
Balmer break signal. The nature of the photoionizing source can 
be inferred from the strength of N iv] and the non-detection of 
N v which also lies in a favourable region of the FORS2 spec- 
trum. Again, by analogy with the arguments already used for the 
Lynx arc (F03; |Binette et al.|2003||Villar-Martm et al.|2004] >, the 
limit on this line ratio can be used to argue for a hot blackbody- 
like source rather than the hard power law spectrum associated 
with an AGN. The absence of detectable radio ( |Kellermann et al.| 
2008 ) or X-ray ( Giacconi et al.|2002] > emission from the source 
at very low flux levels is also evidence against an AGN interpre- 
tation. 

Since our modelling is designed to provide a 'proof of con- 
cept' rather than a detailed SED fit, we use a set of blackbodies 
ranging from 60 to 120 kK as ionization sources. The use of hot 
stellar models having low or zero metallicity would have little 
effect on the gross properties of the predicted spectrum although 
they would be required to fit individual emission line fluxes if 
more observations were available. The effective temperature (the 
tempterature of the black body which produces the same bolo- 
metric luminosity) of such stellar models would be somewhat 



lower (by ss 10-15 kK) than the blackbody temperatures we 
use here. 

The principal requirement of a model is to produce the ob- 
served luminosity of the N iv] doublet. The combination of den- 
sity and ionization parameter, U (defined at the inner edge of the 
nebula) is constrained by the measured upper limit on the source 
size: 



U 



2(H) 



Anr 2 



where Q(H) is the rate of hydrogen ionizing photons radiated 
by the source, r m is the inner radius of the cloud, «(H) is the 
hydrogen number density of the cloud and c is the speed of 
light. Since, for all of our models, the depth of the ionized 
gas is small compared with the observed size of the source, 
we associate r m approximately with the PSF-corrected measure- 
ment from the zgso ACS image of r e (Table [4]) which corre- 
sponds to 580 pc. The observed wavelength of Niv], compared 
with that of the blue-absorbed Lya line, suggests that the dou- 
blet is emitted from a region that is in the high density regime 



(Vanzella et al. ( 2009b i show an enlarge d plot of the spectrum 
of this region) with n e > 10 4 5 cm 4 , see: [Osterbrock & Ferland| 
(2006). We note, however, that in the Lynx arc, the density de- 
rived from Niv] doublet is significantly higher than th e value 
of < 3000 cm~ 3 derived from those of Cm] and Sim] (Villar- 



Martm et al. 2004 



We use the photoionization code CLOUDY version 08.00 
to calculate the models of Hn regions 



(Ferland et al. 



exposed to the blackbody radiation. We assume an ionization 
bounded spherical geometry (i.e., the outer radius of the cloud is 
not defined) which means that the ionizing photon escape frac- 
tion is zero. The nebular gas is assumed to have constant density. 
From the output of the CLOUDY models, which includes both 
nebular continuum and line emission, we compute the theoret- 
ical magnitudes for all the photometric filters to compare with 
the GOODS observations. 

We have calculated a grid of models exploring the param- 
eter space of 7bb of the blackbody, the ionization parameter 
(t/) and the nebular metallicity (Z ne b), investigating how the 
changes in these parameters affect the N iv] luminosity and the 
observed SED. We have assumed the solar relative abundances 
(as defined in CLOUDY) with the realisation that this is un- 
likely to be entirely appropriate at this epoch. In practice, this 
means that our calculated magnitudes will be affected by the 
abundance-dependent strength of lines from individual ions. The 
most important of these is the [Om] 4959, 5007 A doublet that 
is included in the IRAC1 band. Note that, with the range of 
abundance considered here, oxygen is not the dominant coolant 
and so this emission line strength is proportional to abundance. 
Consequently, our principal assumption regarding abundances is 
that of a solar N/O ratio. Guided by the example of the Lynx 
arc, we have performed our modelling assuming that dust is not 
present. However, given the conclusions of the Lye profile mod- 
elling discussed by Van zella et al.| ( [20~09b| l, we do later examine 
the possibility of some dust extinction being present. 

In addition to their fluxes, the emission line equivalent 
widths (EW) are an important constraint because of their de- 
pendence on the strength of the continuum. In our models, the 



' But note that, for the observed line ratio, the N iv] densit y in this 
paper shoul d be 10 s cm -3 , according to the method presented in 
et al. 1 1995| , rather than the quoted 2x 10 4 cirT 3 . 
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continuum under the N iv] line is dominated by hydrogen two- 
photon emission which becomes collisionally de-excited above 
an electron density of ~ 10 4 cm -3 . The observed EW (Table [TJ 
provides an additional constraint on 7\,b, U, n e and Z ne b. 

Guided by the results of x 2 testing of our model grid against 
the observed photometry and N iv] emission line flux, we choose 
a 'fiducial' model that represents our best estimate for the model 
parameters (see Figure [6). Because of our spectroscopic obser- 
vation, high weight is accorded to the luminosity and the equiva- 
lent width of the N iv] doublet. This gives the parameters shown 
in Table 12 
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Fig.5. Observed and modelled SED of GDS18.92-02.7. Top 
panel: the computed nebular spectrum from the fiducial model, 
the Rayleigh-Jeans tail of the ionizing blackbody spectrum and 
the GOODS filter passbands. Bottom panel: the observed and 
computed photometric SED. Additionally VIMOS R (VLT) and 
MIPS 24 pm (Spitzer) magnitudes are included in the SED 
( Santini et al.|2 009). The observed points have l<x error bars or 
are lcr upper limits. The points below Lya have been corrected 
for IGM absorption. Our model does not model the Lya radia- 
tive transfer and so the z'775 point is really an upper limit. Note 
that the Ks band samples almost pure nebular continuum. 



It is significant to note that this is very close to the set of pa- 
rameters as derived for the Lynx arc (F03) based on the fitting 
of many individual emission lines. There are differences in the 
assumed geometry and the particle density but these have only 
minor effects on the emitted spectrum. We have used a higher 
density in order to keep the nebula small enough to be consis- 
tent with the imaging data. The SED derived from this fiducial 
model is shown in Figure [5] and the luminosities of some princi- 
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Fig. 6. As Figure |]but showing the effects of varying the four 
model parameters: upper panel: 7bb, second panel: \og(U), third 
panel: \og(n H ) and lower panel: Z ne b. Note the effect of metal- 
licity on the IRAC1 band (due to [Om]) and the effect of the 
density on IRAC4 due to changes in the Hei 10830A line. To 
improve visibility in this plot, note that the error bars (but not 
the upper limits) are 3<x. 

Table 2. Parameters of the fiducial photoionization model. 



Parameter 


Value 


Unit 




100 


kK 


log(f/) 


-1 




Zneb 


0.05 Z Q 




2(H) 


3 x 10 55 


s- 1 




10 29 


cirr 3 


r\n 


10 21 


cm 



pal emission lines in Table[5] This model produces an H n region 
with an ionization structure and electron temperature distribu- 
tion that results in an radially-averaged electron temperature of 
23 kK. 

In this model, the Ks band samples mostly nebular contin- 
uum while IRAC1 is dominated by very strong [Om] line emis- 
sion resulting in a jump of 1.3 magnitudes between restframe 
wavelengths of 0.33 and 0.54 pm. The IRAC2 band is signifi- 
cantly influenced by Ha and IRAC4 by Hei 10830 A, this latter 
line becoming stronger at higher density (see the third panel of 
Figure 15) . Examination of Figure |6j (which has 3-cr error-bars) 
shows that our fiducial model falls below the low signal/noise 
data from the IRAC3 and IRAC4 filters (see Figure [T). The 
IRAC3 band in our model is dominated by bound-free contin- 
uum and is little affected by changing model parameters within 
the range we have considered. Although the IRAC4 band model 
also falls below the data, the amount can be lessened signifi- 
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Table 3. Fiducial model emission line spectrum. 



Arest 


log(L [erg s -1 ]) 


WW 


Ovi 1035 


40.984 


0.007 


Sim 1207 


41.589 


0.026 


Lyo- 1215 


44.724 


36.011 


Nv 1240 


42.229 


0.115 


Si iv 1397 


42.564 


0.249 


Oiv 1402 


42.416 


0.177 


Nrv] 1486 


42.938 


0.590 


Civ 1549 


44.054 


7.698 


Hen 1640 


43.026 


0.7217 


Cm] 1909 


42.985 


0.657 


[On] 3727 


41.411 


0.018 


Hen 4686 


42.127 


0.0911 


[Om] 5007 


43.936 


5.865 


Ho- 6563 


43.625 


2.867 



Note that no resonance line absorption is accounted for in these 
values. For the doublets, including [O m] 4959, 5007 A, the number 
represents the sum of the two lines. 



cantly by changing the density or adding a higher density com- 
ponent. Given the marginal formal significance of these two dis- 
crepancies we do not consider it justified to attempt to add an 
additional source that would contribute flux at these wavelengths 
without destroying the fit at shorter wavelengths. In order to pro- 
vide enough light from stars, the population would need to be 

~ q2Q09b| > 



old, massive and significantly reddened. IVanzella et al. 



have examined SED fits using both single and fixed stellar popu- 
lation models. These both require an old, > 700 Myr for a single 
and ~ 400 Myr for a mixed, massive, ~ 10 11 M Q , population in 
order to provide the light seen in the IRAC bands. Such an addi- 
tion would completely change the nature of the young, low-mass 
model we are proposing. 

In addition to the requirement to explain the strength of the 
Niv] line, an important constraint is to not over-predict other 
emission lines within the FORS2 spectrum. The absence of N v 
1243 A has already been used to argue against an AGN-like ion- 
izing continuum but it also limits the temperature of the black- 
body to < lOOkK. The Oiv/Siiv complex from 1394-1413 A 
increases in strength at the higher temperatures and is still pre- 
dicted to be significant in our fiducial model (see Table [3]). It 
is not seen in our data but, since we have no a priori informa- 
tion about the abundance of silicon, we do not consider this to 
be sufficient grounds for rejecting the model. The C iv resonance 
doublet is very close to the long wavelength limit of the observed 
spectrum and its detection is dependent on the effective redshift 
which can be influenced, like Lya, by absorption. If it has the 
same redshift as Niv], the observational limit is below the pre- 
diction of our fiducial model by a factor of nearly 20 while if it 
is at the redshift of Lyor it would be undetected (as discussed by 
IVanzella et al.|2009b) . 

3.2. Sample modelling 

The source GDS 18.92-02.7 is the only member of our sam- 
ple that contains a clearly detected second emission line in 
the FORS2 spectroscopy. There is, however, one source with 
groundbased IR spectroscopy that allows the measurement of 
[On], [Om] and H/3 line strengths. This is object CDFS 6664 
(GDS J033233. 33-2 75007 .4) at a redshif t of 3.791 from the 
AMAZE programme ( |Maiolino et al.|2 008 ). 



The presence of significant [On] emission (i.e., comparable 
to H/3) rules out the high temperature and high ionization pa- 
rameter models that we need to explain strong N iv] emission. 
Since all of the spectra reported by Maiolino et al. (2008} do 
show [O n], their use of the locally calibrated R23 parameter and 
other line ratios to estimate metallicity is justified although it 
would not be for a Niv] -emitter, since it was not designed for 
the high ionization parameter cases. While we might expect the 
higher redshift sources to show evidence of hotter ionizing stars, 
the Lynx arc is at a similar redshift to the AMAZE sources but 
shows no [O n] and strong [O m] and N iv] lines. From Figures [2] 
and[3] we see that most of our sample exhibit the rather flat spec- 
tra in f v that is a characteristic of the nebular models. However, 
without the evidence provided by emission line measurements in 
addition to Lya-, it is not possible for us to disentangle the stellar 
and nebular contributions. The relative importance of starlight 
and nebular emission in high redshift starforming galaxies has 
been discussed by Zackrisson et al. (2008 ) who model the effect 
of individual strong emission lines on multi-band SED observa- 
tions. As we have confirmed, both with the Lynx arc and with 
GDS 18.92-02.7 the [Om] lines can carry enough flux to affect 
significantly the magnitudes in even rather broad photometric 
bands and it is in the regime where its strength is proportional 
to oxygen abundance. The presence of a significant excess flux 
in the band that includes these lines can be a hint that we are 
dealing with an SED with a significant nebular component. 



4. Discussion 

The model we employ for GDS 18.92-02.7 is clearly rather rudi- 
mentary and is not expected to produce a detailed match to even 
the somewhat modest data set that we have. Simply summing 
models with a range of model parameters, notably density, would 
be more realistic and would tend to produce a better match to the 
observed SED at the longer wavelengths. However, the general 
properties of this kind of model are clear: the nebular continuum 
is approximately flat in f v and is modulated by the presence of 
strong lines in particular photometric bands. 

Although there appear to be no low redshift analogues of 
GDS 18.92-02.7, the recent discovery of "Green Peas" by the 
Galaxy Zoo project (Cardamone et al. 2009 1 does, however, 
show that there exists a population of compact, strong [Om] 
emitting sources with high specific star formation rates, low 
metallicities and low reddening in the local universe. These 
sources demonstrate that vigorously star forming galaxies can 
exhibit SEDs with very significant nebular contributions to 
broad-band photometric measurements. 

Arguments for thermal rather than non-thermal (AGN) pho- 
toionization were given in section 3.1. It should also be noted 
that the upper limit on the MIPS 24 /zm flux (Table |4] and 
Figure [5]l excludes the steeply-rising SED from warm dust typi- 
cal of AGN of types 1 and 2 ( jMiley & De Breuck|2008] l. Upper 
limits (3<x) are available also on the MIPS 70 and 160 /urn 
fluxes of 2,500 and 33,000 /Jy respectively (Dickinson, M., 
priv. comm.; [Coppin et~aL (2009)) but these values are, unfor- 
tunately, not strong discriminants. The source DLS 1053-0528 
at z = 4.02 discussed by Gli kman et al.| ( (2007] > appears in a list 
of low luminosity type 1 AGN because its emission lines have a 
broad component with a FWHM > 1000 km s _1 . However, aside 
from this characteristic, this object shares emission line and SED 
characteristics more closely with GDS18.92-02.7 and the Lynx 
arc than with typical AGN, notably the great relative strength of 
the nitrogen line with N iv]/ C iv > 1 . 
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The lower limit on the baryonic mass for our fiducial model 
is the sum of the masses of the ionizing stars plus the mass of 
the ionized nebula. To estimate the first of these components we 
use the zero metallicity models of |Schaerer| ( |2002] > which give 
a value of 2(H) per star. We choose his T e g = 85 kK model 
since it produces a similar ionization structure to our fiducial 
model. Such stars have a mass of 50 M G and Q(U) =~ 3 x 10 49 
photons s" 1 . Given the Q(H) required for our model of 3 x 10 55 
photons s _1 , we need a million stars with a total mass of 5 X 
10 7 Mg. The mass of ionized gas from the fiducial CLOUDY 
model is 2.7 X 10 8 M , giving a total minimum mass of 3.2 x 
10 s M Q , more than two orders of magnitude smaller than the 
stellar masses estimated by |Stark et al.| ( |2007] ) and |Wiklind et al.| 
(2008). 

The need to properly include the nebular emission in high 
redshift starforming galaxies is not a new idea, e.g., |Scha erer 
(2002|l; [Zackrisson et al.| l |2008) , but we believe that the Lynx 
arc at redshift 3.4 (F03) and GDS18.92-02.7 at z = 5.6 are the 
only such objects known where the entire observed SED long- 
ward of Lyo; may be contributed by nebular emission. The im- 
plication of this is that most of the stellar flux is emitted in the 
Lyman continuum. As shown by F03, a stellar population with 
a normal IMF that produces sufficient ionizing photons to pro- 
duce the observed emission line fluxes would be easily detected 
in the restframe UV. Unless we have a mechanism, such as a 
dust screen, that selectively absorbs this starlight in our line of 
sight, we are forced to the conclusion that we have a large pop- 
ulation of stars with an effective temperature close to 100 kK. 
The absence of light from cooler stars implies that the cluster 
is either very young or has an intrinsically top-heavy IMF. The 
latter is an expectation for Population III stars. While we have 
no direct method to measure the stellar metallicity in these two 
objects, we do have estimates of the nebular metallicity of m 5% 
solar. This may seem high compared with normal estimates for 



a Pop III environment ( Bromm et al. 2001) but we need to better 



understand the rate of self-pollution within a cluster of a mil- 
lion stars in a deep potential well with temperatures around 80- 
100 kK and lifetimes of just a few million years. 

Our modelling has been carried out without the inclusion of 
dust reddening. It is possible, however, to maintain a reasonable 
fit with the same model parameters other than an increased in- 
trinsic luminosity (and total mass) by including E(B - V) < 0.1 
of Galactic dust for a screen model or E(B- V) < 0.2 for a mixed 
model. The presence of some dust may be necessary to explain 
the Lya profile ( |Vanzella et al.|2009b| ). 

Once the need for a hot, thermal photoionization source is 
identified, it is natural to associate this solely with the presence 
of massive stars. Other sources of thermal radiation could be 
present, e.g., from shocks associated with a high supernova rate 
but, even if the e ntire kinetic luminosity of a core-collapse super- 
nova ~ 10 51 erg (jWoosley & Weaver 1986 1 could be converted to 



as N iv] and He n will also emit strong nebular continuum which 
we show can dominate the Rayleigh-Jeans tail from the ionizing 
stars. The presence of other strong lines falling within a pho- 
tometric band can result in SED features than can be confused 
with jumps such as the Balmer break in conventional stellar pop- 
ulation models. It is important to realise that the detection of a 
luminous emission line such as N iv] necessarily implies a sub- 
stantial nebular contribution to the overall SED and therefore 
effectively excludes models which use purely stellar continuum 
contributions. The large equivalent widths of the emission lines 
can be used to limit the contribution of stars significantly cooler 
than 100 kK to the SED and so constrain the age of the cluster 
or its stellar IMF. 

These highly ionized nebulae appear very different from the 
majority of Hn regions in the local universe due to the higher 
temperature stars and higher ionization parameters. The rest- 
frame UV spectrum is rich with emission lines such as the in- 
tercombination doublets of C, N, O and Si which can be used 
for reliable abundance determinations. 

Determined efforts with large groundbased telescopes 
equipped with efficient red and NIR spectrographs may allow 
us to observe some of the first stellar nucleosynthesis products 
from the young universe. 
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the Niv] line. Daniel Schaerer and Daniel Pequinot gave us much useful mod- 
elling advice. We thank Francesca Matteucci for helpful discussions about the 
chemical evolution of galaxies. AR would like to thank Massimo Stiavelli for 
hosting her visit to STScI to learn and discuss matters pertaining to early star 
formation. In particular, we thank Eros Vanzella for his sterling work on the re- 
duction and analysis of the GOODS/FORS2 spectroscopy and both he and Piero 
Rosati for many inputs to and discussions about this work. 



Lyman continuum photons, this falls far short of the ~ 10 erg 
of ionizing radiation emitted by a massive star during its lifetime 

( |Schaerer|20"0"2] i. 



5. Conclusions 

The possibility of starforming galaxies at high redshift, where 
very hot ionizing stars are expected, having restframe UV 
through MIR SEDs containing significant nebular contributions 
has been discussed for the last decade. The detection of emission 
lines in addition to Lyo- provides, by using photoionization mod- 
els, a way of quantifying this contribution. Compact objects con- 
taining very hot ionizing stars capable of producing lines such 
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